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Abstract. The dependence of the individual and of correlated skeletal variables of Carduelinae and Fringillinae upon bo- 
dy mass (size) was tested by means of multiple regression and single linkage R-cluster analyses. The most obvious re- 
sult is the significantly high degree of body mass dependences of all individual variables within the three functional com- 
plexes feeding, hind limb locomotion and flying. Elimination of body mass as a measure of size decoupled several units 
of correlated variables. Nevertheless the general robustness of units was strong but the removal of "size" had an impact 
on the extent of morphological integration, measurable as agglomeration level of variable units. Body mass (size) in our 
sample thus influenced the quantity rather than the quality of morphological integration. The key feature in fringilline/car- 
dueline skeletal morphology — coevolution of skull and hind limb traits — thus appears to be relatively unconstrained by 


size. 
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1. INTRODUCTION 
1.1. The role of body size in avian morphology 


Body size has been identified as an important factor driv- 
ing the diversification and evolution of organisms. OLSON 
& MILLER (1958, 1999), e. g. quantified size as a strong 
integrative force of morphology, FAIRBAIRN (1997) iden- 
tified its role in sexual dimorphism, KLINGENBERG & 
SPENCE (1997) analyzed its function in life-history evo- 
lution. BONNER (2004) proposed a third size dependent 
rule besides the already well known “weight-strength-rule” 
and “weight-surface-area-rule”. His ^weight-complexity- 
rule" explains that increase and decrease in size of organ- 
isms are strongly related to increase and decrease in com- 
plexity in general, from differentiation of cell types to the 
social organization of organisms. 


These examples implement size as a morphological key 
feature in the organization of life, at all hierarchical lev- 
els from the individual to species (BARBOSA et al. 2000; 
BONNER 2004; CHERNOFF & MAGWENE 1999; PETERS 
1995), 


Size, as the physical magnitude of objects, can be quan- 
tified using either measures or weights. Linear length 
measures are used to seize both the two dimensional area 
as well as the three dimensional volume, whereas weight 
quantifies the proportionality to volume alone. Measures 


widely accepted to represent overall size of organisms in- 
clude also factors like the first principal component 
(BOOKSTEIN et al. 1985; ZELDITCH & FINK 1995; LEISLER 
& WINKLER 2006), landmarks or ratios. BIORKLUND & 
MERILA (1993), e.g. used orthogonal rotation to remove 
size from species mean vectors, LEISLER (1980) calculat- 
ed the cube root of a compound measurement of the body 
core skeleton (sternum length plus pelvis length plus cora- 
coid length) times (sternum width plus pelvis width) times 
height of crista sterni as measure for "body size", follow- 
ing HOERSCHELMANN (1966). 


For the characterization of volume, body mass is the ad- 
equate measure. Several rules explaining the construction- 
al constraints on morphology are related to body mass: The 
"weight-strength-rule" already mentioned above, ex- 
plains why e.g. longer arms, wings and legs have to be 
proportionally thicker than shorter ones, the “weight-sur- 
face-area-rule" describes physiological relations of e.g. gas 
diffusion in large versus small lung surfaces. These are 
geometric growth laws, that explain the change of propor- 
tions — so-called allometric relations — within a biologi- 
cal organism. The underlying process is differential growth 
during development brought about by differences in the 
growth rates of the various parts (WADDINGTON 1966). 


In small bodied flying vertebrates, flight imposes identi- 
cal physical constraints on the flying organism that lead 
to a functional symmetry in body sizes (MAURER et al., 
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2004). The body plan of birds is primarily constrained by 
the weight-strength-rule, that regulates flying ability. In 
extant species 22 kg is the limit for take off; the world- 
wide tiniest bird ıs the Bee Hummingbird, Mellisuga he- 
lenae, found in Cuba and weighing from 1,6-1,9 g (CHAI 
1999). The most successful clade in modern birds, the 
passerines (or perching birds, Passeriformes) have minia- 
turized to an average body mass minimum of 4,5 g in some 
tropical species (Black-faced Flycatcher Warbler, Abro- 
scopus schisticeps, Broad-billed Flycatcher Warbler, 
Tickellia hodgsoni, Cuban Gnatcatcher, Polioptila lem- 
beyei; DUNNING 1993). The heaviest and thus largest song 
bird, the Raven, Corvus corax, from Alaska weighs 1.240 
g. Perching birds successfully evolved and radiated with- 
in any terrestrial habitat —with exception of the polar re- 
gions where their occurrence is physiologically con- 
strained by the weight-surface-area-rule. 


1.2. Key features of the carduelid skeleton 


Our paper focuses on the impact of body mass on the ar- 
chitecture of the avian skeleton, taking the song bird sub- 
families Fringillinae and Carduelinae as a case study. Mor- 
phological diversity in skeletal elements in this phyloge- 
netic entity (clade, monophylum) is mostly confined to 
changes in overall size (BJORKLUND 1991, 1994; VAN DEN 
ELZEN & NEMESCHKAL 1986, 2006; NEMESCHKAL & VAN 
DEN ELZEN 1990, vAN DEN ELZEN & NEMESCHKAL 2006). 
Differentiations in shape are (despite their unquestioned 
ecological importance) only expressed as an insignificant 
portion of the total variation. Whereas three shape axes 
explained only 9,3 % of variance, the size axis account- 
ed for 85,5 % of variance (VAN DEN ELZEN & NEMESCHKAL 
2006). The observed pattern has been identified as a gen- 
eral rule in passerine birds by BJORKLUND (1994, 2006). 
"....making species to be larger or smaller copies of each 
other, all moving up and down a common line of allo- 
metry”. Mapping of skeletal variables along the first prin- 
cipal components (NEMESCHKAL & VAN DEN ELZEN 1990) 
revealed, that morphological diversification in Fringilli- 
nae and Carduelinae focuses on variables of the beak. 
Thus, morphological change in these subfamilies 1s main- 
ly reflected in feeding mechanisms. 


In terms of ecomorphology, a discipline taking the bio- 
logical role of morphological structures under consider- 
ation (BOCK & VON WAHLERT 1965), elements of the avian 
skeleton belong to three different functional complexes: 
skeletal elements of the skull and beak are allied to the 
feeding complex, wing bones, pectoral elements and bones 
of the shoulder girdle represent the flying apparatus and 
skeletal elements of pelvis and legs the hindlimb locomo- 
tion complex. In the carduelid skeleton, nevertheless, units 
of covarying variables (UCCs) contain combinations oth- 


er than proposed by the functional complexes described 
above. These UCCs were detected to reflect units carry- 
ing out common biological roles. Besides a more or less 
undisturbed module of the "flying complex", formed sole- 
ly by wing bones plus elements of the shoulder girdle, the 
variables of the feeding and hindlimb locomotion com- 
plexes were united in one single common cluster, and 
demonstrate the superimposure of feeding on hindlimb lo- 
comotion. 


In search of additional modern approaches in the study of 
morphological integration, CHERNOFF & MAGWENE (1999) 
propose a hierarchical framework for integrative hypothe- 
ses: Morphological variables are judged as integrated at 
the broadest, i. e. the uppermost inclusive level due to their 
covariation with size. At a less inclusive level variables 
covary due to developmental and/or functional associa- 
tions and at still lower levels due to anatomical (spatial) 
associations. Starting out from physics we may assume 
that changes of singular parts from a macroscopical phys- 
ical (e. g., mechanical) system will be ruled by various 
law-dependencies when the whole system changes size 
(e.g., mass). As a consequence, we test in the current study 
the dependence of individual and correlated skeletal vari- 
ables upon body mass taken as a measure of size. 


1.3. Presumptions 


Based upon our knowledge from preceding analyses (VAN 
DEN ELZEN, NEMESCHKAL & CLASSEN 1987, NEMESCHKAL 
& VAN DEN ELZEN 1990, NEMESCHKAL, VAN DEN ELZEN & 
BRIESCHKE 1992, VAN DEN ELZEN & NEMESCHKAL 2006) 
our hypotheses are that: 


l.a significant dependence of skeletal variables on body 
mass exists. lt is expected to differ between variables 
and to be highest 1n beak elements; 

2.thus the three theoretical functional complexes -feed- 
ing, hindlimb locomotion, flying- differ in their size de- 
pendence; 

a. variables of the feeding apparatus — showing 
highest morphological differentiation between 
species — should exhibit the largest amount of 
size dependence, 

b. elements of the hindlimb locomotion complex 
and of the flying apparatus are thought to show 
lesser dependence, as of lesser interspecific 
variability. 

3.units of variables (UCCs) are also expected to be size 
dependent, elimination of “size” is thought to decouple 
variables within units; decoupling is expected to occur 
particularly within “mixed” variable units of beak and 
hindlimb osteometrics. 
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2. MATERIAL AND METHODS 


In 313 specimens, representing 43 species (Table 1), 42 
equivalent skeletal variables (measurements of bone 
lengths, widths and depths) were taken, as described in 
detail in NEMESCHKAL (1999). Twenty measurements stem 
from the feeding apparatus, 13 represent the flying appa- 
ratus and 9 hindlimb locomotion (Table 2). 


Firstly, centroids of 43 species over 42 variables were built 
from the log scaled original data matrix consisting 0f 313 
specimens. This procedure was essential, because body 
mass data were available as species means only. Morpho- 
logical variation between species was then quantified by 
the total variance of species centroids. Body mass varia- 
tion in the actual fringilline-cardueline sample ranges from 
a maximum of 54 g in the Hawfinch Coccothraustes coc- 
cothraustes to a minimum of 9 g in neotropical siskins 
(Table 1). To test hypotheses about size impact, linear re- 
gression analysis was applied. Log transformed species 
means of body mass were taken as predictor variable and 
each of the 42 skeletal measurements as criterion variable. 
The 42 resulting residuals are used as the variables under 
study — the variables corrected for body mass (= variables 
of which size was partialized out; BCM). The coefficients 
of determination between body mass and original variables 
(Table 2) were tested for significance using random boot- 
strap (1000 replicates each; for computer programme 


coefficient of determination 


1 3 5 7 ot Ws ze = 





package see NEMESCHKAL (1999)). They are figured as 
profiles of variance (Figs 1, 2). Single linkage R-cluster 
analyses were chosen to figure correlations of original da- 
ta (Figure 3- BO; based upon the variable intercorrelation 
matrix between species means for original data) and cor- 
relations data corrected for body mass (Fig. 3- BCM; vari- 
able intercorrelation matrix between species means for 
body mass adjusted data). Units of correlated variables 
were taken into account, when they agglomerated at the 
uppermost quartile (r? > 0.924) in the original dataset and 
reappeared in the size corrected dataset again. E-units are 
stable clusters with topological identity of variable posi- 
tions in both the R-cluster analyses of original data and 
the body mass corrected data set. S-units are less stable 
units built of identical variables but with topological dif- 
ferences between analyses of original and size corrected 
variables. 


3. RESULTS 


The most obvious findings are on the one hand the sig- 
nificantly high degree of body mass correlation of all in- 
dividual variables (Fig. 1) and on the other hand the ro- 
bustness of several covarying variable units found in both 
the original and the "residual" data set, adjusted for body 
mass (Fig. 2). 


21 23 
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skeletal characters 


Fig. 1. Dependence of 42 skeletal characters on body mass as shown by coefficients of determination. 


28 Renate VAN DEN ELZEN & Hans L. NEMESCHKAL: Impact of Body Mass on Morphological Integration in Avian Skeletons 


variances and residuals 


VN  ——A- 
Ud NJ AGLI 





{ 3 5 7 9 1 13 15 17 19 21 23 25 27 29 31 3393553: Booms 
skeletal characters 


— cB) == BCM 


Fig. 2. Comparison of variance profiles. Variance of original characters (BO) occur in the upper line and residuals for characters 
corrected for body mass (BCM) at the bottom. 
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Figs 3. Single linkage R-cluster analysis of original data (BO; left) and of body mass corrected data (BCM; right). Black squares 
indicate measures of the pectoral girdle and forelimbs, black triangles mark skull measurements and open circles measures of the 
pelvic girdle and hindlimbs. Clusters are explained in the text. 
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3.1. Influence of body mass on individual variables 


All 42 variables are correlated with body mass at a high 
significance level (p< 0.001) as revealed by random boot- 
strap procedures. The degree of correlation with “size”, 
measured by coefficients of determination (COD) between 
(log transformed) original data of species means and body 
mass, is highest in coracoid length (variable number 35; 
COD=0,91), post orbital length (12; 0,90) and equal at 
COD= 0,89 in tibiotarsus width (27), total skull length (16) 
and synsacrum length (21). The largest amount of body 
mass independence exhibit carpometacarpus (41) and ul- 
na lengths (39; COD- 0,34 and 0,47) as well as tar- 
sometatarsus length (No 28) and width (29; COD- 0,70, 
resp. 0,67) and proximal end width of the humerus (38; 
COD= 0,71). Skull variables in general reveal the high- 
est correlation with body mass, hindlimb and forelimb 
variables indicate a larger variation in their correlation 
with body mass. 


A comparison of variance profiles of original and size cor- 
rected skeletal variables (Fig. 2) confirms and. strength- 
ens these findings: The relatively higher variances in the 
original variables appear smoothened in the size correct- 
ed residuals (BCM; bottom line in fig. 2), indicating that 
the largest amount of total variance between species is due 
to a “size factor", measurable by body mass. This espe- 
cially holds true for skull variables (skeletal variables 
1—20) and less for hindlimb and forelimb measures (skele- 
tal variables 24—29, 37—42). 


3.2. Influence of body mass on united variables 


The single linkage cluster analysis of original data (Fig. 
3) extracted two main units fused at a high level, several 
disintegrated variables and one isolated variable set (El). 
The module BOC 1 is composed of the variables of the 
flying apparatus, mainly its “engine-supporting” part, in- 
cluding measures of sternum, coracoid and scapula. Wing 
measures that represent "flight performance" are highly 
disintegrated, only humerus length (37) and carpometacar- 
pus width (42) are included. The second main module 
BOC 2 combines skull and hindlimb measures. 


The two main modules BOC 1 and BOC 2 are composed 
of two different types of subsets: Robustness is observed 
in nine smaller units (E1—E9), where only two to three 
variables are coupled. Similar structure is maintained af- 
ter correction for size in four units (S1—S4), including three 
to 17 variables as well as including several of the stable 
units (E4—E9). Whereas, for example, femur length (24) 
and tibiotarsus length (26) are next neighbours in the S4 
cluster in BO, in BCM tibiotarus length (26) and tar- 
sometatarsus length (28) are agglomerated at the highest 


level. The larger S1 and S2 units also include some of the 
smaller, topologically stable E units: E4-E6 always oc- 
cur within S1, E8-E9 within S3, S4 is always part of SI. 


BOCs identified in the R-cluster analysis of original da- 
ta are to a great extent also found in the analysis of body 
mass corrected data (Fig. 3). Again in the body mass cor- 
rected data analysis, a main cluster composed of engine- 
supporting variables from the flight complex (S2 in BOCI, 
fig. 3) opposes a main cluster structured by a mixture of 
skeletal elements of the feeding and hind limb locomo- 
tion apparatus characters. The main differences to the orig- 
inal data analysis are the rather low agglomerative level 
(coefficients of determination in the "size-corrected" 
dataset being two to three times lower than in the origi- 
nal dataset), reduction of the involved variable units and 
some changes in variable topology. BOC2, that was com- 
posed of S1 (including S4), S3, E2, E3 and E7 in the orig- 
inal dataset, comprises in the size corrected data set only 
S1 (including S4) and E3, the remaining components — 
units S3, E2 and E7 — are clustered at lower agglomera- 
tion levels. 


Subcluster S1 differs in its hierarchical structure of ag- 
glomeration levels of E units and arrangement of variables 
within S clusters. In S2 sternum length (30), keel length 
(32) and keel depth (31) are differently arranged, in S3 
post orbital length (12), and within S4 femur-, tibiotarsus 
and tarsometatarsus lengths (24, 26, 28) changed their po- 
sition along branches. 


In both illustrations of single linkage R-cluster analyses 
(Fig. 3), EI shows up as a more or less disintegrated unit. 
It consists of measures of ulna (39) and carpometacarpus 
lengths (41), which are tied together at a relatively high 
agglomeration level, especially in the analysis of "size- 
corrected" data, in which both measures are combined 
with carpometacarpus width (42). 


4. DISCUSSION 
4.1. Size and morphological integration 


Do these results confirm or contradict our predictions and 
expectations of hypotheses on morphological integration? 
Corroborating our first assumption, a general body mass 
dependence was identified at a high significance level for 
all individual characters within all three classical function- 
al complexes. “Size” could thus be confirmed as a strong 
integrative force in cardueline and fringilline skeletal mor- 
phology. Variability was, as assumed in our second hy- 
pothesis, relying to a greater extent on body mass in skull 
and beak measurements than in characters allied to the oth- 
er functional complexes. Elimination of “size”, decoupled, 
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Table 1. List of species under study, their body masses and information sources. Sources: DU= Dunning 1993; RE= 
labels of specimens at ZFMK, Bonn; RO= McLean 1988; WA= Ward 2003. 








Species body mass in g sources 
Fringilla montifringilla 24 DU 
Fringilla coelebs 21,4 DU 
Coccothraustes coccothraustes 54 DU 
Pyrrhula pyrrhula 21,8 DU 
Pinicola enucleator 56.4 DU 
Loxia curvirostra 40,38 DU 
Chloris chloris 215 DU 
Chloris sinica 31.3 DU 
Acanthis flammea I5 DU 
Carduelis cardnelis 15,6 DU 
Linaria cannabina 15,3 DU 
Agriospiza flavirostris 15.4 DU 
Rhodospiza obsoleta 295 DU 
Spinus spinus 14,5 DU 
Spinus cucullatus 9 RE 
Spinus barbatus 16,6 DU 
Spinus xanthogaster 127 DU 
Spinus yarrellii 95 RE 
Spinus psaltria 2.5 DU 
Serinus pusillus 11,6 DU 
Serinus serinus Lb DU 
Serinus syriacus 12.1 BE 
Serinus canicollis 13,8 DU 
Serinus citrinella 12 DU 
Alario alario 11,8 RO 
Serinus canaria 12:97 DU, RE 
Ochrospiza reichenowi IES RE 
Ochrospiza atrogularis 11.4 DU 
Ochrospiza leucopygia 9 RE 
Ochrospiza mozambica 10,6 DU 
Ochrospiza dorsostriata 14,4 DU 
Ochrospiza xanthopygia 10 RE 
Crithagra gularis 15,5 DU 
Crithagra striolata 22.4 DU 
Crithagra albogularis 259 DU 
Crithagra donaldsoni buchanani 24 DU 
Crithagra sulphurata 19:7 DU 
Serinops flaviventris 16,3 DU 
Crithagra mennelli [5:3 RO 
Crithagra burtoni 29:9 DU 
Dendrospiza hyposticta 1 DU 
Dendrospiza scotops 15,4 RO 


Pseudochloroptila totta 13,4 WA 


Bonner zoologische Beiträge 56 (2007) 31 


as predicted in the third assumption, several units of cor- 
related variables. Nevertheless, contradicting our expec- 
tations, displacement of variables from units or units from 
larger modules was low and robustness of units strong. 
Body mass (size) itself in our sample thus influenced the 
quantity rather than the quality of morphological integra- 
tion and did not affect the key feature in fringilline/car- 
dueline skeletal morphology — coevolution of skull and 
hindlimb variables. 


According to classical hypotheses on morphological in- 
tegration, "size" 1s the uppermost integrative level, at deep- 
er agglomerative levels, traits covary due to developmen- 
tal or functional associations and at still deeper levels due 
to spatial adjacency (CHERNOFF & MAGWENE 1999) re- 
spectively morphological neighborhood (ALPATOV & 
BOSCHKO-STEPANENKO 1928). Indeed, in our study, re- 
moval of the first hierarchical layer "size", mostly effect- 
ed the degree of morphological integration, measurable as 
agglomeration level of character units. Structure and re- 
lation of character compositions however remained to a 
great extent constant and robust. For variable units at the 
body mass reduced level, three types of trait correlations 
are identifiable: S-modules are subunits agglomerated at 
higher integrative levels within BOC] and BOC2, the orig- 
inal, not body mass reduced, modules. E-modules that are 
not composed of characters from immediate morpholog- 
ical neighborhood or partial overlap (E4, E5, E9) can be 
assigned to two categories: Classical functional units (E2, 
E3, E6 and E7) are with one exception (E7) fused at low- 
er agglomerative levels (below 0.6 COD) than the units 
of biological roles El and ES (above 0.8 COD). Our find- 
ings thus corroborate the integrative hypotheses quoted 
above, that anatomical association is a primary integra- 
tive level in morphology. 


4.2. Interpretation of the observed character mod- 
ules 


EBLE (2005) classifies four different kinds of modules in 
morphology: 


1. modules due to structural relations, 

2. modules of pleiotropic genotype-phenotype mappings, 

3. developmental units (gene expressions, domains of epi- 
genetic dynamics, regions with localized allometric 
growth) and 

4. functional units. 


We would like to extend the definition of the last catego- 
ry by modules accomplishing common biological roles 
(NEMESCHKAL et al. 1992) and thus differentiate between 


4.a. mechanical functioning (like kinetics) and 
4.b. biological functioning (like morphological traits used 
in €. g. feeding). 


Units E4 (postorbital width and internasal width, 9, 10), 
E5 (skull width and interorbital maximal width, 6, 7) and 
E9 (premaxilla length and total skull length, 11, 16) may 
be explained by their immediate morphological neighbour- 
hood and /or partial overlap (“rule of neighbourhood", AL- 
PATOV & BOSCHKO-STEPANENKO 1928; structural relations 
EBLE 2005). Most of the other results allow for mechan- 
ical functional explanations, few for explanations of their 
biological roles. E2 is composed of two skull characters 
(caudal length of the mandible and quadratojugal length, 
3, 13) and a pelvis measure (synsacrum distal length, 23). 
Whereas the first two skull variables are functional coun- 
terparts and act in seed husking, their —also body mass in- 
dependent- correlation with a synsacrum length measure 
comes unawares. E3 combines two mandible depths (4, 
5) with quadratum length (19), characters functioning in 
the kinetism of the avian skull, especially in lateral bill 
movements enhancing seed husking (NUUENS & ZWEERS 
1997). Also units E6 (skull height and premaxilla depth, 
15, 20) and E7 (femur and tibiotarsus width, 25, 27) are 
easier to explain by constructional constraints than bio- 
logical roles. E6 depicts the interdependence of skull and 
beak height, E7 that of leg diameters, in ES.Only units E1 
and ES where measures of the flying apparatus (car- 
pometacarpus (41) and ulna lengths (39)) respectively bill 
tip are tied together (length of pars symphysialis (1) with 
dentary length of the mandible (2)) allow for an interpre- 
tation of biological roles: The first unit comprises the 
feather-carrying bones of the wing, responsible for ma- 
noevrability in flight per se. Unit ES represents "the bill 
tip grasp’, a feeding tool well developed in carduelid 
finches (NEMESCHKAL et al. 1992). 


S4 with the lengths of femur (24), tibiotarsus (26) and tar- 
sometatarsus (28) obviously reflects a character pattern of 
concerting limb dimensions as guided by developmental 
processes during ontogeny (NEMESCHKAL 1999). 


Comparison of interspecific correlation patterns 
(macroevolutionary pattern, this study Fig. 3) with infra- 
specific analyses of variable correlations (microevolution- 
ary level, NEMESCHKAL & VAN DEN ELZEN 1994, fig. 7) re- 
veal that six modules (E3, E4, E5, E7, E8 and E9 of this 
study) correspond to units also found within carduelid 
finches and three modules (E5, E7 and E9) even to units 
occurring within the phylogenetically distant pigeons. At 
the microevolutionary level NEMESCHKAL et al. (1992) pro- 
vided evidence for clade specificity of variable units (mod- 
ules) and their correspondence to expressions of develop- 
mental control genes (NEMESCHKAL 1999). Morphologi- 
cal modularity might generally be seen as product of de- 
velopmental modularity, because morphological patterns 
of organization emerge in ontogeny (EBLE 2005). Conse- 
quently, each of the extracted variable complexes might 
additionally also correspond to developmental or body 
plan modules. 
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Table 2. Skeletal variables under study and their coefficients of determination. BO represent original data, BCM resi- 
duals for the body mass reduced variables. 


Skeletal variables 





Skull measures 


Coeff. of determination 


Ol mandible, pars symphysialis length 0.7014 
02 mandible, dentary length 0.7641 
03 mandible, pars caudalis length 0,7642 
04 mandible, dentary depth 0.6924 
05 mandible, pars caudalis depth 0,7883 
06 skull width 0,8239 
07 interorbital, maximal width 0,8353 
08 interorbital, minimal width 0,6704 
09 postorbital width 0.8048 
10 internasal width 0,7856 
11 premaxilla length 0.8240 
12 postorbital length 0,86 
13 quadratojugal length 0,7955 
14 orbital width 0,7888 
15 skull height 0.7280 
16 skull length 0,8610 
17 palatinum width 0,6176 
18 pterygoid length 0,0932 
19 processus orbitalis quadrati length 0,7988 
20 premaxilla depth 0,7911 
Pelvic girdle and hindlimbs 

21 synsacrum length 0,7748 
22 synsacrum width 0,7076 
23 synsacrum distal length 0.7464 
24 femur length 0.8355 
25 femur width 0,7781 
26 tibiotarsus length 0,7675 
27 tibiotarsus width 0.7562 
28 tarsometatarsus length 0,6180 
29 tarsometatarsus width 0.3034 
Pectoral girdle and forelimbs 

30 sternum length 0,8413 
31 keel depth 0,6897 
32 keel length 0,7632 
33 sternum, distal width 0,6617 
34 sternum depth 0.7395 
35 coracoid length 0,8401 
36 scapula length 0,8114 
37 humerus length 0,7644 
38 humerus, proximal end width 0,6246 
39 ulna length 0.3253 
40 ulna wıdth 0,6235 
4] carpometacarpus length 0.2535 
42 carpometacarpus width 0,6454 
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How genes and developmental pathways co-operate is best 
documented for the avian beak. It is constructed by mul- 
tiple facial prominences (HELMS et al. 2005, HELMS & 
BRUGMANN 2007), where the frontal nasal mass (FNM), 
lateral nasal prominences (LNP), and maxillary promi- 
nences (MXP) comprise the upper beak and the mandibu- 
lar prominence (MDP) forms the lower beak. The identi- 
ties of facial prominences are specified early in the neu- 
ral crest stage and may involve homeobox genes such as 
Hox (COULY et al. 1998) or MSX (BROWN et al. 1997, Wu 
et al. 2006). 


In metazoa, at least 17 signal transduction pathways op- 
erate to activate or repress different genes at distinct times 
and places in the embryo. Five predominate in early em- 
bryonic development: the Wnt, TGF-beta, Hedgehog, 
RTK, and Notch pathways. Five more are used in late de- 
velopment, and seven more in the functions of differen- 
tiated cells (GERHART 1999). 


Diversity in beak shapes of Galapagos finches, e. g., is 
known to originate from the activity of a morphogenetic 
bone protein (BMP4 of the TGF-beta family; ABZHANOV 
et al. 2004) and calmodulin, a protein that binds and trans- 
ports calcium ıons (CaM; ABZHANOV et al. 2006). Where- 
as BMP4 stimulated growth of beaks along two dimen- 
sions, — it produced deeper and wider beaks and thus ex- 
plained the linkage in the variation of these traits, — the 
authors found CaM as regulator of beak length. Thus in 
Galapagos finches, beak length develops independently 
from width and depth due to two different factors that lead 
to changes in growth along different dimensions 
( ABZHANOV et al. 2006). In ducks, chickens and cockatiels, 
BMPA is also differently expressed (WU et al. 2004, Wu 
et al. 2006). At late stages of development, chicken and 
duck embryos had two localized growth zones in the 
FNM, which melted in the chicken, but stayed separated 
in the duck. Ducks, moreover exhibited a wider FNM and 
more activity in another growth factor (fibroblast growth 
factor 8). In cockatiels a thicker FNM increased in a dif- 
ferent direction and the mandibular prominence (MDP) 
was suppressed. BMP4 was involved in all species in me- 
diating activity in all localized growth zones. Experimen- 
tal overexpression in BMP4 altered beak shapes among 
all species and beak curvation was induced by asymmet- 
ric growth activity in a facial prominence. BMP4 was al- 
so found to be responsible for regulating a homeobox 
gene: As it increased the expression of MSXI, the authors 
suggested that MSXI activity is regulated by the BMP 
pathway. 


4.3. Conclusions 


Under the assumption that BMP4 is a main regulator for 
the expression of developmental control genes in the avian 


skeleton, this protein, for instance, can be seen as a key 
growth factor accounting for size and shape variation in 
our skeletal variables. A developmental growth factor can 
also explain the strong size dependence of single variables, 
but cannot elucidate the interspecific correlation patterns 
of variable units. Thus in the light of current knowledge 
on gene expression variability and research in develop- 
mental genetics the classical view on morphological in- 
tegration and genesis of modularity might be too one- 
sided. 


WAGNER et al. (2007) offer a broader approach to modu- 
larity. Summarizing empirical evidence, they discern be- 
tween two levels: The elements the modules consist of and 
their kinds of interaction. Elements vary from molecular 
level (nucleotides) to morphological traits and their con- 
nections from physical to dynamical and statistical. They 
define three kinds of modules: Variational, functional and 
developmental modules, but stress a duality, as modules 
may belong to several categories at the same time. Units 
observed in the present study on cardueline/fringilline 
skeletal variables fall in the category of variational mod- 
ules, as their measure is statistical. For the interpretation 
of variational modules the authors offer several answers: 
identical developmental origin, similarity of covariation 
and gene expression territories, and pleiotropic effects. So 
our variable units (or character complexes; NEMESCHKAL 
et al. 1992) can be defined not only as a duality but a plu- 
rality: according to their operability they are functional 
modules, according to their mode of exploration they are 
variational modules, and according to their genesis devel- 
opmental modules. 


The open problem is, whether modules arise through the 
action of natural selection or because of biased mutation- 
al mechanisms (WAGNER 2007). Both neutral models and 
models based on natural selection are offered: e. g. LYNCH 
(2007) favors the neutral model and stated that “... emer- 
gent biological features such as complexity, modularity, 
and evolvability, all of which are current targets of con- 
siderable speculation, may be nothing more than indirect 
by-products of processes operating at lower levels of or- 
ganization." We do not follow the author in refusing nat- 
ural selection and external evolutionary forces, but as- 
sume, that selection via fitness promotes certain pheno- 
types derived from their genetic basis and processes in on- 
togeny. Under the viewpoint of developmental constraints 
as internal evolutionary forces, we like to generalize the 
findings of Wu et al. (2004, 2006) on the origin of beak 
shape in birds. The authors offer a more reconciling con- 
clusion, that morphological diversity may be achieved by 
modulating prototypical molecular and cellular modules. 
Transforming growth factors may mediate the range, lev- 
el, and duration of locally enhanced growth, thus provid- 
ing a spectrum of morphological designs for selection (Wu 
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et al. 2006). Also WAGNER et al. (2007) arrive at the con- 
clusion that mutational processes (internal processes) fa- 
vor the origin of modularity and selection pressures (ex- 
ternal processes) reinforce the mutational bias. 


Acknowledgements. We thank H. Brieschke and H. Classen for 
collecting data. S. Rick, S. Neulinger and S. Klein patiently dig- 
itized data and prepared figures. Internal institutional discussions 
and the comments of two anonymous referees greatly improved 
the manuscript. 


REFERENCES 


ABZHANOV, A., KUO, W. P., HARTMANN, C., GRANT, B. R., GRANT, 
P. & TABIN, C. J. 2006. The calmodulin pathway and evolu- 
tion of elongated beak morphology in Darwin's Finches. Na- 
ture 44: 563—567. 

ABZHANOV, A., PROTAS, M., GRANT, B. R., GRANT P. & TABIN, 
C. J. 2004. BMP4 and morphological variation of Beaks in 
Darwin's Finches. Science 305: 1462-1465. 

ALPATOV, W. W. & BOSCHKO-STEPANENKO, A. M. 1928. Varia- 
tion and correlation in serially situated organs in insects, fish- 
es and birds. American Naturalist 62: 409-429. 

BARBOSA, A., BARLUENGA, M. & MORENO, E. 2000. Effects of 
body mass on the foraging behaviour of subordinate Coal tits 
(Parus ater). The Ibis 142: 428-434. 

BJORKLUND, M. & MERILA, J. 1993. Morphological differentia- 
tion in Carduelis finches: adaptive vs. constraint models. Jour- 
nal of Evolutionary Biology 6: 359—373. 

BJORKLUND, M. 1991. Patterns of morphological variation among 
Cardueline finches (Fringillidae: Carduelinae). Biological 
Journal of the Linnean Society 43: 239—248. 

BJORKLUND, M. 1994. Process generating macroevolutionary pat- 
terns of morphological variation in birds: a simulation study. 
Journal of Evolutionary Biology 7: 727—742. 

BJORKLUND, M. 2006. Trait evolution, morphological integration 
and functional complexes. Acta zoologica Sinica 52 (Supple- 
ment): 494—497. 

BOCK, W. & von WAHLERT, G. 1965. Adaptation and the form- 
function complex. Evolution 19: 269—299, 

BONNER, J. T. 2004. Perspective: the size-complexity rule. Evo- 
lution 58: 1883-1890. 

BOOKSTEIN, F. L, CHERNOFF, B., ELDER, R., HUMPHRIES, J., SMITH, 
G. & STRAUSS, R. 1985. Morphometrics in evolutionary bi- 
ology. Academy of Natural Sciences, Philadelphia. Special 
Publication 15. 

BROWN, J. M., ROBERTSON, K. E., WEDDEN, S. E. & TICKLE, C. 
1997. Alterations in Msx1 and Msx2 expression correlate with 
inhibition of outgrowth of chick facial primordia induced by 
retinoic acid. Anatomy and Embryology 195: 203—207. 

CHAI, P. 1999, Species account Melisuga helenae. P. 671 in: DEL 
Hoyo, J., ELLIOTT, A., SARGATAL, J. (eds.) Handbook of the 
birds of the world. Volume 5. Barn-owls to Hummingbirds. 
Lynx Edicions, Barcelona. 

CHERNOFF, B. & MAGWENE, P. M. 1999. Morphological integra- 
tion: forty years later. Pp. 319—353 in: OLSON, E.C., & MILLER, 
R.L. (eds.) Morphological integration. University of Chica- 
go Press. 

CouLY, G., GRAPIN-BOTTON, A. COLTEY, P., RUHIN, B. & LE 
DOUARIN, N. M. 1998. Determination of the identity of the de- 
rivates of the cephalic neural crest: incompatibility between 


Hox gene expression and lower jaw development. Develop- 
ment 125: 3445—3459. 

DUNNING, J. B. 1993. CRC Handbook of Avian Body Masses. 
CRC Press, Boca Raton. 

EBLE, G. J. 2005. Morphological Modularity and Macroevolu- 
tion: Conceptual and Empirical Aspects. Pp. 221—238 in: 
CALLEBAUT W, RASSKIN-GUTMAN, D. (eds.) Modularity. Un- 
derstanding the Development and Evolution of Natural Com- 
plex Systems. MIT Press, Cambridge, Massachusetts, London, 
England. 

FAIRBAIRN, D. J. 1997. Allometry for sexual size dimorphism: 
Pattern and process in the coevolution of body size in males 
and females. Annual Review of Ecology, Evolution, and Sys- 
tematics 28: 659—687. 

GERHART, J. 1999. 1998 Warkany lecture: signaling pathways in 
development. Teratology 60: 226—239. 

HELMS, J. A. & BRUGMANN, S. A. 2007. The origin of species- 
specific facial morphology: the proof is in the pigeon. Inte- 
grative and Comparative Biology 47: 338—342. 

HELMS, J. A., CORDERO, D. & TAPADIA, M. D. 2005. New insights 
into craniofacial morphogenesis. Development 132: 851-861. 

HOERSCHELMANN, H. 1966. [Allometric investigations on rump 
and wing of snipes (Charadriidae und Scolopacidae)]. Allo- 
metrische Untersuchungen an Rumpf und Flügel von Schnep- 
fenvögeln (Charadriidae und Scolopacidae). Zeitschrift für 
zoologische Systematik und Evolutionsforschung 4: 209-317. 

KLINGENBERG, C. B. & SPENCE, J. R. 1997. On the role of body 
size for life-history evolution. Ecological Entomology 22: 
55—68. 

LEISLER, B. 1980. Morphological aspects of ecological special- 
ization in bird genera. Ecology of birds 2: 199-220. 

LEISLER, B. & WINKLER, H. 2006. Interaction between morpho- 
logical integration and ecological factors. Acta Zoologica Si- 
nica 52 (Supplement): 498-501. 

Lyncu, M. 2007. The fragility of adaptive hypotheses for the ori- 
gins of organismal complexity. Proceedings of the National 
Academy of Sciences 104: 8597—8604. 

MAURER, B. A., BROWN, J. H., DAYAN, T., ENQUIST, B. J., ERNEST, 
S. K. M., HADLY, E. A., HASKELL, J. P.. JABLONSKI, D., JONES, 
K.E., KAUFMAN, D.M., Lyons, S.K., NIKLAS, K.J., PORTER, 
W.P., Roy, K., SMITH, F.A., TIFFNEY, B. & WiLLIG, M. R. 2004. 
Similarities in body size distribution of small-bodied flying 
vertebrates. Evolutionary Ecology 6: 783—797. 

MCLEAN, G. L. 1988. Roberts’ Birds of Southern Africa. New 
Holland Publishers, London 

NEMESCHKAL, H. L. & ELZEN, R. VAN DEN 1990. [Functional 
complexes in carduelid skeletons — a morphometric approach]. 
Funktionskreise im Skelettsystem der Cardueliden. — Ein mor- 
phometrischer Ansatz. Pp. 37—41 in: ELZEN, R. VAN DEN (ed.) 
Current Topics in Avian Biology. Proceedings of the Interna- 
tional Centennial DO-G meeting, Bonn. 

NEMESCHKAL, H.L., ELZEN, R. VAN DEN & BRIESCHKE, H. 1992. 
The morphometric extraction of character complexes accom- 
plishing common biological roles: avian skeletons as a case 
study. Journal of Zoological Systematics and Evolutionary Re- 
search 30: 201-219. 

NEMESCHKAL, H. L. & ELZEN, R. VAN DEN 1994. Macrospecific 
character coupling — an analysis of morphological integrati- 
on in avian skeletons. Zoologische Jahrbücher: Abteilung für 
Systematik, Ökologie und Geographie der Tiere 121: 265-302. 

NEMESCHKAL, H. L. 1999, Morphometric correlation patterns of 
adult birds (Fringillidae: Passeriformes and Columbiformes) 
mirror the expression of developmental control genes. Evo- 
lution 53: 899—918. 


Bonner zoologische Beiträge 56 (2007) 35 


Nuyens, F. W. & ZWEERS, G. A. 1997. Characters discriminat- 
ing two seed husking mechanisms in finches (Fringillidae: 
Carduelinae) and estrildids (Passeridae: Estrildidae). Journal 
of Morphology 232: 1-33. 

OLSON, E. C. & MiLLER, R. L. 1958. Morphological integration. 
University of Chicago Press. 

OLSON, E. C. & MILLER, R. L. 1999. Morphological integration. 
University of Chicago Press. 21d edition. 

PETERS, R. H. 1993. The ecological implications of body size. 
Cambridge studies in ecology. Cambridge University Press. 

VAN DEN ELZEN, R. & NEMESCHKAL, H. L. 2006. Variabilität im 
Skelettbau der Cardueliden — ökologisch, phylogenetisch, ad- 
aptiv oder konstruktionsmorphologisch dominiert? Ökologie 
der Vógel 25: 231-247. l 

VAN DEN ELZEN, R., NEMESCHKAL, H. L. & CLAREN, H. 1987. 
Morphological variation of skeletal characters in the bird fa- 
mily Carduelidae: I. General size and shape patterns in Afri- 
can canaries shown by principal component analyses. Bon- 
ner Zoologische Beiträge 38: 221-237. 

WADDINGTON, C. H. 1966. The Modular Principle and Biologi- 
cal Form. Pp. 20-35 in: Kepes G. (ed.) Module, Proportion, 
Symmetry, Rhythm. Vision and Value series. George Braziller, 
New York. 

WAGNER G. P, MEZEY J & CALABRETTA, R. 2005. Natural selec- 
tion and the origin of modules. Pp. 33—49 in: CALLEBAUT W. 
& RASSKIN-GUTMANN, D. (eds.) Modularity. Understanding 
the development and evolution of natural complex systems. 
MIT Press, Cambridge, Massachusetts, London, England. 

WAGNER, G. P. 1996. Homologues, natural kinds and the evo- 
lution of modularity. American Zoologist 36: 36—43. 

WAGNER, G. P., PAVLICEV, M. & CHEVERUD, J. M. 2007. The road 
to modularity. Nature Reviews Genetics 8: 921—931. 

WAGNER, G. P. & ALTENBERG, L. 1996. Complex adaptations and 
the evolution of evolvability. Evolution 50: 967—976. 


WARD, V. L. 2003. Safring results. Pseudochloroptila totta. South 
African Bird Ringing Unit Report 855. Avian Demography 
Unit, Cape Town. 

Wu, P. , JIANG, T.-X., SIKSAWEANG, S., WIDELITZ, R. B. & 
CHUONG, C.-M. 2004. Molecular shaping of the bill. Science 
305: 1465-1466. 

Wu, P., JANG, T.-X., SHAN, J.-Y., WIDELITZ, R. B. & CHUONG, 
C.-M. 2006. Morphoregulation of avian beaks: comparative 
mapping of growth zone activities and morphological evolu- 
tion. Developmental Dynamics 235: 1400-1412. 

ZELDITCH, M. L. & Fink, W. L. 1995: Allometry and develop- 
mental integration of body growth in a piranha, Pygocentrus 
nattereri (Teleostei: Ostariophysi). Journal of Morphology 
223: 341-355. 


Authors’ addresses: Renate VAN DEN ELZEN (correspon- 
ding author), Section Ornithology, Zoologisches 
Forschungsmuseum Alexander Koenig, Adenauerallee 
160, D-53113 Bonn, Germany; E-Mail: r.elzen.zfmk@uni- 
bonn.de; Hans L. NEMESCHKAL, Department of Theoret- 
ical Biology, University of Vienna, Althanstrasse 14, A- 
1090 Wien, Austria; E-Mail: Hans.Leo.Nemeschkal@uni- 
vie.ac.at 


Received: 21.02.2008 

Revised: 25. 06. 2008 

Accepted: 09. 09. 2008 
Corresponding editor: M. Schmitt 


ZOBODAT - www.zobodat.at 


Zoologisch-Botanische Datenbank/Zoological-Botanical Database 

Digitale Literatur/Digital Literature 

Zeitschrift/Journal: Bonn Zoological Bulletin - fruher Bonner Zoologische Beitrage. 
Jahr/Year: 2009 

Band/Volume: 56 


Autor(en)/Author(s): Elzen Renate van den, Nemeschkal Hans. L. 


Artikel/Article: The Impact of Body Mass on Morphological Integration in Avian Skeletons 
(Aves, Fringillidae; Carduelinae, Fringillinae) 25-35 


